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INTRODUCTION 


On  approach  to  ah  airless  pla:i«.*»  such  as  the  moon*  the 
slowdown  from  approach  speed  to  a  aoft-landing.tfpeed  must  be 
accompUohed  entirely  by.an  on-board  propulsion  system*  since 
no  aerodynamic  drag  Is  available  fsr  this  purpose.  The  minimum 
approach*  or  unslowed  impact  sf^ed  is  the  escape  velocity*  which 
for  the  mcon  is  7780  feet  per  second.  An  additional  speed  of  the 
order  of  1000  feet  per  second  is  usually  desirable  in-order  to 
.horten  the  transit  time  and  to’tailor that'tirhe-so-thc  lunar  landing 
j  occurs  when  the  moon  is  visible  from  chosen  locations  on  the  earth. 
By  a  soft  landing  is  meant  a  speed  of  about  10  to  15  feet  per  second* 
impl^dng  a  maximum  deceleration  on  impact  of  10  to  15  earth  g. 

\  I 

Thus,  the  retro  engine  system  must^tak©  out  some  9000  feet  per 
^  '  second  of  velo'cUy*  ,  The  Ji^^ortancc  of  efficient  design  of  tho  descent 
maneuver  and  of  the  ret^o  system  can  be  emphasized  by  stating  that 
approximately  75  percent  of  the  injected  spacecraft  weight  is  devoted 
to  the  retro  system.  Of  this  75  percent*  90  percent  is  fuel.  A 
small  fractional  change  in  retro  system  weight  therefore  represents 
a  very  large  fractional  change  in  payload  weight. 
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The  desigQ  and  optimixation  of  the  descent  trajectory 
mlalmum  retro  system  weight  is  considered  here  fo^.an 
vehicle.  The  design  procedure  is  not  a  theoretical  one  with^O|;S£^ 
example,  'n' propulsion  stages  and  a  calculus  of  variations^^timl- 
zation.  Instead,  the  major  characteristics  of  the  retro  ays^xr^r^ 
chosen  from  a  practical  hardware  consideration  of  simple,  rbrlihole 
engines  and  Instruments  which  exist  today.  Once  chosen,  th^^p^^ 
tical  limitations  of  these  components  impose  constraints  on  1 
syctem,  and  it  is  within  these  engineering  constraints  that^ 
optlmizaticn  procedure  is  performed. 

In  Part  Z  of  this  paper,  the  physics  of  the  procedure  is 
emphasized,  not  the  mathematics.  In  fact.  Part  I  illustrates  the 
procedure  in  general  under  the  simplifying  assumptions  that 
1)  the  trajectories  are'ver  tical,  without  velocity  and  attitude 
direction  being  considered  as  variables,  and  2}  the  gravity  field 
is  constant.  Under  these  assumptions,  the  mathematics  are  very 
sir.ple  oir.ee  the  equations  of  motion  are  integrable  in  closed  farm 
and  a  linear  de.*ivative  technique  performs  the  error  analysis.  The 
effects  resulting  from  '  onsidering  trajectory  and  thrust  attitude 
angles  off  vertical  as  variables  are  mentioned  briefly  In  Parti* 

The  detailed  method  of  handling  the  optimization  procedure  for 
machine  computation  In  this  more  complicated  case  appears  in 
Part  Jl. 
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PART  I.  ENGINEERING  DESIGN  AND  OPTIMIZATION 
OF  TRAJECTORIES 


1.  MAJOR  DESIGN  CONSIDERATIONS 

It  Is  well  known  that  for  ideal  efficiency  a  slowdown  In  a 
gravity  field  should  be  made  at  maximum  thrust  for  minimum  time 
and  that  the  maneuver  should  start  at  an  altitude  where  essentially 
zero  velocity  is  reached  at  precisely  the  moon's  surface.  Such  a 
maneuver  requires  resistance  to  unlimited  acceleration  and  perfect 
sensing  and  control  of  the  trajectory*  In  a  practical  attempt  to 
approach  this  maneuver«  a  high-level,  controllable  thrust  engine 
is  indicated.  This  cohtrollabiUty. demands  a  liquid  propellant 
engine.  For  comparable  engine  weight,  howe'^er,  present  liquid 
propellant  engines  have  one -third  the  thrust  level  of  solid  propellant 
engines.  This  lower  thrust  level  is  less  efficient  in  itself  and  also 
requires  ignition  at  a  much  higher  altitude,  therefore  necessitating 
much  heavier  instrumentation.  For  an  unmanned  vehicle  at  dis¬ 
tances  close  to  the  moon,  the  greater  accuracy  and  faster  action 
of  a  control  system  which  responds  automatically  to  local  sensing 
devices  is  preferable  to  one  relying  on  tracking  information  and 
command  from  the  earth. 

For  the  major  portion  of  the  slowdown,  a  high  constant- 
thrust,  8olid*propellant  engine  is  desirable.  Since  this  engine  is 
not  controllable,  altitude  and  velocity  dispersions  will  result  at 
burnout.  The  nominal  burnout,  therefore,  must  occur  at  an  altitude 
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high  enough  to  cneure  that  a  second*  controllablei  llquid^propellant 
engine*  responding  to  local  sensing*  can  remove  the  remaining* 
relatively  small  velocity  before  impact  occurs* 

Figure  1  illustrates  typical  descent  trajectory  behavior*  In 
this  figure  and  the  following  fljures*  however*  the  trajectory  is^not 
pictured  in  space*  since*  although  the  ordinate  is  altitude*  the  abscissa 
is  velocity*  downward  to  the  right*  upward  to  the  left.  The  actual 
trajectories  are  to  be  considered  vertical*  ^or  off-vertical  approaches* 
the  ordinate  can  be  regarded  as  range  and  the  abscissa  as  range  rats* 
This  approximation  is  good*  since  in  the  case  of  off-vertical  trajec¬ 
tories  with  thrust  direction  closely  opposing  velocity  direction*  the 
angle  of  approach  changes  very  jittle  until  very  low  speeds  are 
reached. 

In  Figure  1*  the  dashed  line  to  the  right  represents  an  approach 
trajectory  without  slowdown*  the  velocity  increasing  slighUy  as  the 
vehicle  falls  to  the  moon.  The  solid  curve  represents  the  trajectory 
during  main  retro  slowdown.  As  the  trajectory  proceeds*  there  is 
less  altitude  loss  for  two  reasons:  1)  the  average  speed  is  lower* 
and'2}  for  a  near-constant  thrust  engine*  the  acceleration  is  greater 
since  mass  is  decreasing.  The  dashed  line  in  the  center  represents 
a  fall  from  a  midvalue  of  velocity.  Note  that  now  the  increase  of 
speed  for  a  given  altitude  drop  is  greater  than  for  the  unslowed 
vehicle.  Since  the  average  speed  is  less*  more  time  is  spent  in 
falling  through  a  given  height*  and  more  time  in  a  gravity  field 
means  more  speed  increase.  This  increase  in  speed  must  always 
be  taken  out  if  a  soft  landing  is  to  be  accomplished. 
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Toward  the*  left  in  Figure  1*  the  fall  starting  at  a  low  down¬ 
ward  velocity  exhibit#  a  large  speed  increase  in  a  smaU  altitude  drop. 
Indicated  at  the  extreme  left  is  the  absurdity  of  having  such  a  large 
impulse  engine  that  the  velocity  at  burnout  is  upward.  The  vehicle 
coasts  upwardi  slowing  down«  then  falls  from  a  greater  height  acquir¬ 
ing  even  more  speed  before  impact. 

One  possible  simple  concept  for  the  second  vernier  engine 
•maneuver  is  illustrated  by  the  vertical  dot-dash  curve.  Here  the 
second  engine  simply  holds  a  thrust  level  equal  to  lunar  weight  so 
that  the  descent  from  main  retro  burnout  is  madc^t  the  constant  low 
landing  speed.  Such  a  maneuver  is  good  and  will  be  adopted  at  the 
very  end  of  the  descent,  the  final  drop  of  a  few  tens  of  feet.  But 
the  maneuver  is  poor  for  the  whole  descent  from  main  retro  burnout 
for  two  reasons:  1)  the  velocity  dispersion  at  main  burnout,  as  will 
be  shown,  is  too  large  to  accept  as  a  landing  speed  range,  and  2}  the 
altitude  dispersion  is  large  enough  that  even  if  the  burnout  speed  were 
corrected  by  the  vernier  engine  to  the  low  landing  speed  soon  after 
main  burnout,  the  fall  at  this  low  constant  speed  through  the  remain¬ 
ing  altitude  at  vernier  thrust  equal  to  lunar  weight  would  be  very 
'Costly  in  vernier  engine  fuel. 

The  desired  trajectory  after  main  retro  burnout  results  from 
the  same  considerations  which  governed  the  main  retro  slowdown. 

It  is  best  to  drop  as  quickly  as  possible  to  an  altitude  and  velocity 
condition,  sensed  by  the  Instruments  on  board.  Operation  of  the 
vernier  engine  at  its  highest  thrust  level  will  then  bring  the  vehicle 
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to  essentially  zero  velocity  at  zero  altitu^le*  Such  a  trajectory 
would  not  use  the  vernier  engine  during  the  first  phase  of  the  descent 
following  main  retro  burnout.  It  is  desirable^  however#  to  use  the 
vernier  engine  for  attitude  control  during  the  entire  powered  descent# 
since  stabilization  by  spinning  is  awkward  for  the  .sensing  instruments# 
and'Since  the  possible  misalignment  of  the  main  retro  thrust  is  suffi¬ 
cient  to  require  a  sizable  correcting  moment  to  ensure  proper  attitude 
behavior. 

Intermittent#  om^off#  operation  of  the  vernier  engine  presents 
a  problem  in  that  continuous  attitude  control  is  needed#  and  also 
presents  a  greater  reliability  problem  than  does  continuous  operation 
which  controls  thrust  level  through  a  range  between  minimum  and 
maximum  thrust.  Is.  attitude*  control  and  reliability#  then#  the 
decision  is  to  fall  after  main  retro  burnout  at  a  mlnimxxm  vernier 
thrust  level  of  somewhat  less  than  lunar  weight.  The  known  weight 
of  the  vehicle  together  with  the  maximum  to  minimum  thrust  ratio  of 
a  practical#  liquid  propellant  vernier  engine  thus  establishes  the 
available  thrust  levels  of  the  vernier  engine  in  pounds. 

The  higher  the  maximum  to  minimum  thrust  ratio  obtainable# 
the  more  efficient  is  the  vernier  descent.  Beyond  a  certain  value 
of  this  ratio#  however,  a  variable  injector  area  Is  necessary#  whereas 
the  more  limited  ratio  obtainable  by  pressure  control  alone  has  the 
advantage  in  reliability  of  a  proved  design. 

Thus,  from  consideration  of  the  major  characteristics  of 
engines  and  Instruments,  a  two-englnc  retro  system  la  selected. 
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An  is  shown  in  Figure  2,  the  first  xcain  solid  propellant  engine 
removes  . the  bulk  of  the  velocity  at  a  high  constant  level  of  thrust. 

A  second*  controlVible*  three>chamber»^  liquid  propellant  vernier 
engine  controls  vehicle  attitude  during  the  main  retro  engine*  then 
operates  at  its  minimum  thrust  following  main  retro  burnout  \mtil  a 
sensed  velocity~altitude  comblxxation  is  reached  such  that  maximum 
vernier  thrust  ^operation  brings  the  vehicle  to  a  speed  of  perhaps 
10  feet  per  second  at  a  few  feet  above  the  lunar  surface.  Then 
operation  at  a  thrust  level  equal  to  lunar  weight  holds  this  low  speed 
until  impact*  or  allows  the  vehicle  to  drop  in  free  fall  for  the  last 
few  feet.  Thus  the  major  characteristics  of  the  retro  system  are 
decided. 

For* instrumentation*  this  engine  system  requires  a  range 
sensor  for  triggering  ihe  main  retro  at  the  proper  altitude  and  a 
precision  range  sensor  and  velocity  magnitude  and  direction  sensor 
for  monitoring  the  vernier  descent  phase  both  in  thrust  program  and 
in  attitude  slaving  of  the  thrust  vector  to  the  velocity  direction. 

Following  a  quantitative  evaluation  of  the  main  retro  errors 
and  dispersions  to  be  expected*  it  now  is  a  matter  of  choosing  the 
exact  aivision  of  energy  removal  between  the  first  and  second  engine 
such  that  the  total  retro  system  weight  is  least.  The  optimization 
procedure  will  answer  the  question*  '*At  what  nominal  speed  whould 
the  burnout  of  the  main  retro  engine  occur?" 


2.  MAIN  RETRO  ENGINE  DISPERSIONS 


Tor  a  vertical  trajectory*  In  whlcli  the  thruet  line  is  coniidered 
perfectly  aligned  to  the  vertlc^  velocity^  there  are  five  eburcee  of 
deviations  which  produce  dispersions  Jn  the  burnout  alutude  and 
velocity^  two  sources  of  Initial  conditions  and  three  sources  of  engine 
conditions.  Under  the  assumptions  of  a  constant  gravity‘field»  Table  1 
gives  the  equations  for'^1)  the  main  retro  descent*  2)  the  dispersions 
incurred  during  this  descent*  and  3)' the  minimum  thrust  and  maxi¬ 
mum  thrust  phases  of  the  vernier  descent  considering  vehicle  mass 
constant  during  the  vernier  descent. 

Figure '3  illustrates  these  error  contributions*  and  the  result¬ 
ing  ellipse  represents  dispersions  in  altitude  and  velocity.  An  alti¬ 
tude  marking  Inst^rumcnt  provide  ii  the  signal  for  ignition  of  the  main 
retro  engine^  Error  In  tMs  measurement  results  In  error  of  bum- 
•out  altitude  only*  in  a  constanyr^ravlty  field;  a  closely  approximate 
truth.  Because  of  imperfect  injection  or  midcourse  correctlpn 
conditions*  an  error  in-izxitial  velocity  is  also  present.  This  error 
produces  a  dispersion  at  burnout  in  both  altitude  and  velocity*  a 
correlated  dispersion  which  contributes  to  the  inclination  of  'the 
dispei'sion  ellipse.  Note  here  that  errors  In  velocity  detected  at 
midcourse  of  the  same  order  of  magnitude  as  the  velocity  dispersions 
of  the  main,  retro  engine  are  better  left  uncorrected  in  the  midcourse 
phase*  since  the  particular  velocity  error  to  be  contributed  by  ^e 
main  retro  engine *l8  as  yet  unknown.  Darger  velocity  corrections 
desired  may  be  efficiently  corrected  at  midcourse.  It  is  pertinent 
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to  mention  here  that  the  vernier  engine  logically  performs  tlie  short 
period  i^dcourse  maneuver  to  correct  injection  errors  which  would 
result  priinarily  in  error  in  the  lunar  landing  location* 

A  variation  In  the  amount  of  main  retro  propellant  loaded  and 
of  sliver  behavior  at  burnout  determines  the  exact  amount  of  propellant 
burned  and  ^erefore  contributes  to  variations  in  the  burnout  velocity* 
but  has  a  net'Usibie  effect  on  the  burnout  altitude*  Variation  in'the 
specific  impulse*  a  variation  in  chemistry*  affects  total  impulse  by 
affecting  thruot  level*  Such  a  variation  also  produces  changes  in 
both  burnout  velocity  and  altitude*  Finally*  variations  in  burning 
rate*  as  caused  by  grain  shape  deviations  and  more  importantly  by 
varying  temperature  of  the  propellant  as  affected  by  ambient  launch 
pad  and  transit  conditions  cause  changes  in  thrust  level  with  essen¬ 
tially  no  change  in  total  impulse.  Their  effect  is  therefore  essen¬ 
tially,  limited  to  burnout  altitude  with  little  offect  on  burnout  velocity* 
Although  not  shown  in.Figure  3*  if  there  is  a  variation  in 
total 'VeMcle  weight  at  main  retro  ignition  due  to  variation  In  the 
midcourse  correction  fuel  consumed*  the  effect  is,  similar  to  that 
caused  by:  variation  in  both  main  engine  fuel  burned  and  in  thrust 
level*  This  dispersion  is  handled  automatically*  however*  If  the 
nominal  ignition  weight  is  taken  as-corresponding  to  previous  com¬ 
plete  consumption  of  the  midcourse  fuel  allowance*  The  heavier 
weight  requiring  more  vernier  descent  fuel  Is  due  to  having  more 
vernier  fuel  available* 
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Table  1«  Vertical  Descent  Equations 


Definition  of  Symbols  (In  order  of  use) 

(  }  Value  at  start  of  main  retro 
'  *o 

(  Value  at  end  of  main  vr'^TO,  start  of  minimum  thrust 

vernier  phase 

(  Value  at  end  of  minimtun  thrust  vernier*  start  of  maximum 

thrust  vernier  phase 

V  Velocity*  feet  per  second*  negative  when  downward 

h  Altitude*  feet 

c  Exhaust  velocity  tr  G)*  feet  per  second 

"Wp  Propellant  weight  consiimed,  pounds 

ini.4al  spacecraft  weight*  po\-.r.ds 
g  Lunar  surface  gravity  acceleration,  54  3Zft/sec^* 

consld*ired  constant 

G  Earth  surface  gravity  acceleration*  32.Zft/sec^ 

T  Thrust*  pounds 

b  Ratio  of  maximum  vernier  thrust  to  lunar  spacecraft  weight 

a  Ratio  of  minimum  vernier  thrust  to  lunar  spacecraft  weight 

t  Time*  seconds 

m  Spacecraft  mass*  slug<i(«  considered  constant  during  vernier 

descent 

Main  Retro  Descent 

r  (w„/w  ;,g/G)' 

Vj  =  Vo  -  cl  log(l  -  Wj,/Wo)  +  - 
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Table  1 ;  (continued) 


*‘1 '  h 


{a-Wp/w^-; 


*+V<=- 


p/W^-)  log(l  -  w^/wj 
,  «W/WJ(g/G)11 


‘i'tiuin  'R'itro  Burnout  Dl«per«iont 


i(T/W„) 

avj  -  4V^  +  - -  .  — ,  ° 


+  eW  /W^ 
P  o 


cV„  VJ-W 


■*■  “GT/^ — - 


Vernier  Minimxim  Thruat 
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Table  1  (conUnued) 

Vernier  Maximum  Thruat  Peacent  (ending  at  rero  velocity  and  ^ltltude} 
2g{b  -  Db^ 

=-V2(b/b-l) 

■I  “  /max  T  ^ 


Combined  Vernier. Minimum  and  Maximum  Thrv>5t  Petc>*nt 


3.  VERNIER  rUEL  CONSUMPTION 

For  the  purpose  of  illustrating  the  optlmixation  procedure, 
j  Figure  d'shcws  the  latter  part  of  the  descent  trajectory  under  the 
assumption  of  a  given  lunar  impact  speed,  at  a  larger  scale  than 
in  the  previous  figures.  The  upper  solid  curve  shows  the  end.  of 
the  nominal  main  retro  descent,  and  the  ellipse  shoves  the  dispersions 
about  the  nominal.  The  lower  solid  curve  represents  the  descent 
trajectory  using  n«axim\im  vernier  thrust.  The  vehicle  must  never 
be  allowed  to  descend  into  the  altitude-velocity  region,  below  and 
to  the  right  of  this  curve,  since  insthat  region  the.maximxun  thrust 
of  the  vernier  engine  is  too  low  to  slow  the  vehicle  to  a  soft  landing 
speed  before  impact,  regardless  of  the  amount  of  vernier  fuel 
available. 

Loci  of  constant  vernier  fuel  consumption  are  shown  by  each 
dashed  curve.  These  curves  indicate  altitude-velocity  combinations 
from  which  vehicle  descent  ui\der  the  prescribed  minimum  and  then 
maximum^vemler  thrust  phases  consumes  a  given  constant  amount 
of  fuel.  As  indicated  in  Figure  4,  the  higher  constant  fuel  curves 
represent  more  vernier  fuel  consumed.  From  the  nominal  burnout 
condition,  the  vehicle  drops  along  the  dot-dash  curve  at  minimum 
vernier  thrust  to  the  solid  curve.  At  this  point,  altitude  and 
velocity  sensors  command  maximum  vernier  thrust,  and  the 
vehicle  follows  this  .naximxim  thrust  trajectory  to  approximately 
zero  velocity  at  zero  altitude.  There  is  a  point  on  the'boundary  of 
the  dispersion  ellipse,  toward  the  upper  right  of  the  ellipse  shown 
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in  Figure  4,  from  which  a  deacent  consumes  the  maximum  amount  of 
vernier  fuel  corresponding  to  the  particular  nominal  burnout  condition 
Illustrated. 

4.  OPTIMIZATION  PROCEDURE 

Consider  the  effect  of  moving  the  nominal  burnout  point  and 
therefore  the  dispersion  ellipse  to  the  left  from  a  location  near  the 
line  of  peaks  of  the, constant  vernier  fuel  curves.  This  move  repre-^ 
sents  &  main  retro  burnout  at  a  slower  speed  and  therefore  requires 
more  main  retro  fuel.  From  the  shapes  of  the  constant  vernier  fuel 
cxirves,  it  is  apparent  that  this  ch:uige  also  requires  more  vernier 
fuel.  Clearly*  such  a  main  burnout  region  is  undesirable.  Consider* 
however*  Increasing  the  main  burnout  velocity*  moving  the  ellipse  to 

I 

the  right  of  the  peaks  of  the  constant  vernier  fuel  curves.  Again  more 
vernier  fuel  is  required*  but  now*  with  less  velocity  decrement  during 
the  main  retro  descent*  less  main  retro  fuel  is  needed.  Thus  in  this 
region  of  main  burnout  speed*  a  compromise  exists*  indicating  the 
possibility  of  an  optimum  main  burnout  velocity  such  that  the  total 
descent  fuel  required  is  least. 

Figure  5  shows  the  change  in  main  retro  fuel*  in  vernier  fuel* 
and  in  total  fuel  as  a  function  of  nominal  main  burnout  velocity. 

Again  the  lunar  approach*  or  rather  impact  speed*  Is  considered 
constant.  It  is  clear  front  Figure  5  that  an  optimum  main  burnout 
speed  does  in  fact  exist*  one  for  which  vernier  fuel  is  at  a  minimum* 
more  importantly  one  for  which  total  descent  fuel  is  at  a  minimum. 
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To  avoid  misunderstanding*  It  Is  emphasized  that  Figure  5  does  not' 
represent  absolute  weight  values.  The  zeros  o£  the  different  curves 
are  different;  only  changes  are  shown. 

Of  course*  it  is  least  total  retro  system  weight  thatrls  desired* 
rather  than  Just  descent  fuel  weight.  As  nominal  main  burnout  velocity 
changes*  so  does  main  retro  ignition  altitude*  and  the  initial  altitude 
of  vernier  descent.  These  altitude  changes  have  their  effect  on  the 
range  capability*  and  therefore  weight*  of  the  sensing  instrximents. 
More  importantly*  dry  engine  weight  changes  along  with  fuel  capacity. 
Taking  these  changes  Into  account*  the  curves  of  Figure  5  may  be 
considered  to  represent  change  in  fuel-plus-dry-engine  weight  and 
change  in  total  retro  system  weight*  including  sensing  Instrumenta¬ 
tion.  Regarding  Figure  5  in  this  light*  it  demonstrates  that  an 
optimxxm  division  of  velocity  removal  between  the  first  and  second 
engine  has  been  found  such  that  total  retro  system  weight  is  least* 

The  optimization  has  been  accomplished;  the  retro  system  is  designed. 

5.  UVUNCH  FL.EX1BIL.ITY 

With  the  retro  system  defined*  the  velocity  decrement  of  the 
main  retro  engine  and  the  fuel  capacity  of  the  vernier  engine  are 
fixed.  For  an  approximately  constant  transit  time  to  the  moon* 
determined  by  the  locations  on  the  earth  of  the  launch  site  and  the 
tracking  facility*  the  lunar  approach  speed  will  vary*  being  in 
general  greatest  at  lunar  apogee  and  least  at  perigee.  In  addition* 
the  varying  lunar  declination  causes  transit-time*  and  therefore 
-approach-speed  variations  for  earth-fixed  visibility  requirements. 
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For  flexibility  of  launch  time  throughout  the  montli*  the  vernier  fuel 
capacity  io  choaen  to  handle  the  highest  approach  speedy  and  the 
optimization  is  performed  for. this  moat  severe  condition*  For  condi¬ 
tions  other  than  apogee>  the  nominal  main  burnout  will  occur  at  a 
speeo  lower  than  that  shown  In  Figure  5  corresponding  to  least  total 
retro  system  weight.  As  nominal  main  burnout  speed  decreases^ 
Figure  5  shows  that  vernier  fuel  required  at  first  decreases,  then 
rises  again.  The  descent  is  still  successful  until  such  a  low  burnout 
speed  is  reached  that  the  vernier  fuel  reqxilred  is  again  equal  to  that 
for  the  design  apogee  condition.  Within  this  range  of' nominal  burnout 
speeds  and  therefore  of  approach  speeds#  which  vary  during  the  month, 
the  veliicle  retro  system  is  flexible. 

Should  greater  extension  of  the  range  at  the  low  speed  end  b6 
desired,  It  can  bo  gained  at  some  cost  in  simplicity.  Depending  on 
injection  error,  some  vernier  fuel  may  be  available  for  an  efficient 
midcourse  velocity  correction.  Alternatively,  since  near  perigee  the 
injection  velocity  is  lower,  more  weight  can  be  injected,  and  a 
greater  ignition  weight  would  raise  the  nominal  burnout  speed.  The 
selection  of  a  main  retro  engine  off-loaded  by  a  few  pounds  would 
accomplish  the  same  result. 

6.  EFFECT  OF  OFF-VERTICAL  TRAJECTORIES  AND  ATTITUDES 

It  has  been  stated  that  the  altitude  and  downward  velocity  of  the 
vertical  trajectory  may  be  effectively  considered  as  range  and  range 
rate  for  off-vertical  trajectories.  This  is  especially  true  if  the  thrust 
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vector  U  always  closely  aligned  In  opposition  to  the  velocity-  vector. 
In  that  case>  the  trajectory  direction  changes  very  slowly  during 
most  of  the  descent  and  rotates  to  vertical  only  in  the  last*  very 
low  velocity  phase.  Moreover*  as  is  assumed  here*  if  the  main 
retro  thrust  direction  is  held  constant  at  the  initial  velocity  direc¬ 
tion  at  ignition*  the  trajectory  changes  appreciably  from  its  Initial 
direction  only  In  the  las'  few  seconds  of  main  retro  operation. 
Velocity  sensing  and  therefore  servo  aligning  of  the  thrust  vector 
to  the  instantaneous  vslocity  vector  during  the  main  retro  descent 
is  difficult  because  of  the  dense*  highly  ionised  exhaust  products 
of  the  solid  engine.  The  far  lower  thrdst  level  and  therefore  less 
dense  exhaust  Jet  of  the  vernier  engines  is  not  expected  to  be 
troublesome  to  sensorsi  and  velocity  information  during  the  vernier 
phase  to  servo  the  thrust  in  opposition  to  the  velocity  direction  is 
relied  upon.  Should  the  velocity  direction  be  more  than  perhaps 
45^  off-vertical  at  the  end  of  main  retro*  however*  the  sensing 
of  this  ^ct  and  s^se^uent  slaving  of  the  vehicle  (that  Is*  thrust) 
attitude  to  this  direction  could  result  in  such  an  extreme  attitude 
that  the  velocity  sensors  no  longer  can  perform  their  service  for 
lack  of  view  of  the  moon's  surface. 

To  avoid  this  possibility*  the  full  problem  must  be  examined 
in  which  the  velocity  and'thrust  attitude  direction  are  additional 
variables.  In  this  case*  the  equations  of  motion  are  no  longer 
integrabla  in  closed  form  and  machine  computation  becomes 
Inevitable.  The  planar  curves  shown  for  example  in  Figure  4 


in  this  ditcuf  tlon  of  vortlc&l  tr&Joctorioi  become  three-dimeneional 
sur^ceA;  the  dispersion  ellipse  becomes  a  dispersion  fooUbU.  The 
optimisation  proceduro  is  very  similar  to  that  already  described,  but 
the  mathematical  techniques  are  more  complicated  and  warrant  a 
treatment  of  ^eir  own  In  Fart  XI  of  this  paper*  Suffice  it  to  say  here 
that  the  effect  of  attitude  errors  and  off-vertical  trajectories  show 
importance  largely  on  the  low  end  nominal  burnout  velocity  cutoff. 
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PART  U.  A  DIGITAL  COMPUtER  PROGRAM  FOR 
OFF-VERTICAL  TRAJECTORIES 


1.  MAIN  RETRO  DISPERSION  ELLIPSOID 

As  in  the  case  ol  the  vertical  trajectories  discussed  in 
Part  Ii  the  Urst  step  in  computation  for  the  off-ve^tical  case  is 
an  evaluation  of  the  dispersions  at  main  retro  burnout.  For  ver¬ 
tical  trajectoriesi  there  were  five  sources  of  burnout  dispersions, 
two  in  initial  conditions^and  three  in  engine  behavior.  To  these 
must  now  be  added  three  new  dispersion  sources.  The  first  is  an 
error  in  the  direction  of  the  initial  velocity  at  main  retro  ignition. 
The  second  is  an  error  in  the  vehicle,  and  therefore  thrust,  atti¬ 
tude  in  that  it  is  not  perfectly  aligned  in  opposition  to  the  initial 
velocity.  The  third  error  results  from  the  type  of  thrust  attitude 
control  assumed  here. 

In  Part,  1  no  local  sensing  of  the  instantaneous  velocity 
vector  during  main  retro  operation  was  assumed  because  of  the 
difficulty  of  sensing  through  the  dense  exhaust  Jet.  In  Part  U 
the  thrust  attUuCo  iu  similarly  considered  constant  at  the  initial 
velocity  direction  during  the  major  portion  of  main  retro  opera¬ 
tion.  In  order  to  set  up  the  computations  for  the  full  range  of 
po  -  «ble  variables,  thus  broadening  the  scope  of  the  program,  it 
ie  now  assumed  that  at  a  few  seconds  before  main  burnout,  local 
sensing  and  consequent  slaving  of  the  thrust  attitude  to  the  instant¬ 
aneous  velocity  vector  occurs.  The  third  additional  error,  then, 
is  a  bias  error  in  this  continuous  alignment  of  thrust  to  the  chang¬ 
ing  velocity  vector. 

The  eight  sources  of  dispersion  are  listed  as  follows,  in 
the  order  in  which  they  are  used  in  the  following  analysis: 

Error  in  thrust  attitude  from  initial  velocity  vector 
Error  in  thrust  attitude  from  instantaneous  chang¬ 
ing  velocity  vector 

Variation  in  velocity  magnitude  at  ignition 
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At  Variation  in  thrust  due  to  propellant  temperature 
variation 

Variation  in  specific  impulse  (therefore  also  thrust) 
At  I  Variation  in  burning  time*  equivalent  to  variation 

-in  propellant  burned 
Ah^  Variation  in  ignition  altitude 
Ay^  Variation  in  velocity  direction  at  ignition 

Bvcause  of  these  eight  perturbations*  there  are  main  retro 
dispersions  in  alticude*  h*  velocity  magnitude*  V*  and  now  also  in 
velocity  direction*  y*  thus  forming  c.  three-dimensional  space. 

£aeh  of  these  three  dispersions  is  a  variable  dependent  on  the 
eight  perturbations  which  are  considered  as  random  variables. 

The  random  variables  are  assumed  to  be  normally  distri¬ 
buted  and*  pairwise*  mutually  independent.  Nominal  trajectories 
eve  run  on  the  IBM  7090  for  the  various  combinations  of  initial 
velocity  orientation  (relative  to  the  horisontal)  and  retro  terminal 
velocity  magnitude:  i.e.  *  10  or  20  cases*  as  required.  For  each 
of  these  nominal  conditions*  the  variables  as  defined  are  perturbed 
one  at  a  time*  and  the  corresponding  deviation  from  the  uomlnal 
is  obtained.  Thus*  the  set  of  independent  varlabler  is  reduced  to 
three  mutually  dependent  variables  at  the  end  of  the  first  retro 
phase.  These  three  are  assumed  to  be  Jointly  normally  distributed. 
A  confidence  ellipsoid  is  then  defined  by  the  resulting  three- 
dimensional  normal  distribution.  In  order  to  make  this  a  99  per¬ 
cent  confidence  Interval*  the  quadratic  form  exponent  is  rotated 
to  its  normal  form  and  the  chi-squared  distribution  for  three  de¬ 
grees  of  freedom  is  employed  to  establish  the  constant  value  for 
the  exponent.  This  approach  also  establishes  the  eemi-major  axes 
and  the  orientation  of  the  ellipsoid. 

The  equations  of  motion  for  V  >  V^,  where  the  thrust 
vector  is  fixed  in  space  are  given  below  and  the  geometry  Is  shown 
in  Figure  6. 
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These  equations  of  motion  are  employed  until  V  a  V  r 
which  time  the  following  equations  are  employed  (see  Figure  7). 
The  initial  conditions  for  these  equations  are  the  values  of  y, 
V  .  V  ,  etc.  at  V  =  V_. 


Figure  7 
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(16)  y 
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This  8et'o£  equations  is  solved  simultaneously  until  V  s 
--  some  prescribed  main  retro  terminal  velocity.  The  time  to 
reach  Vj  is  then  observed,  and  the  perturbation  trajectories  are 
run  .  3ing..this  yalu>j<o£ 

From  the  above  information  and  under  the  hypothesis  of 
small  perturbations,  the  variablei^  h,  V,  and  defined  as  linear 

combinations  of  the  eight  perturbniion  variables  prevtoUsly  listed. 
The  perturbation  notation  is  dropped  for  convcnlence--'i.e. , 
etc. 

The  various  coefficier.ts  are  the  respective  terminal,  per 
unit  perturbations  from  the  iceal  trajecvory,  a^ ,  and  d^'. 

(18)  h  =  ajQo  +  =  +  agti  +  *8^  o 

(19)  V  =  b,e^  b^Gj  +  bjT  +  b^c  +  bgtj  +  b^V^  +  b^ho  +  bgY„ 

(20)  Y  =  djG^  +  djGj  +  djT  +  d^c  +  dgt,  +  dgV^  +  d7h^  +  dgY  ^ 

Since  each  of  the  independent  variables  of  equations  (18), 
(19)>  and  (20)  is^Caussiau  wUh  known  parameters,  the  parameters 
of  the  joint-distribution  can  be  determined  for  h,  V,  and  y*  Thus, 
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(22) 


i=  1 


is  1 


where  x.  =  6  *  ©j,  T  .  .  •  etc. 


=  atandard  deviation  of  the  random  variable 


Phv  =  E 


PVv  =  E  •'A 

is  1 


PYh  =  S  Vi‘"*f/'"v"'h 


s  correlation  coefficient  between  the 
*^nin  ^ 

m  and  variable 
Tn,n  =  b,  V#  Y 

In  terms  of  equations  (22)  through  (26)  and  under  the  assumption 
that  the  joint  distribution  is  Caussian,  this  probability  density  is 
defined  as 


(27)  p(h,V,Y)  s 


mW'  |A|V»  - 


along  with  the  covariance  matrix  equation  (28). 


7V  = 


PVh  "-v 

‘■vh'v  "h 


PhV  ®’h  ®'v 

‘‘yV'y  'v 


“  determinant  of  the  covariance  matrix 

=  cofactor  of  the  m*^  row  and  the  n***  column  term 
=  h.  V,  V 

This  distribution  is  considered  in  terms  of  0  meane  for  convenience. 

Once  the  confidence  ellipsoid  is  established^  ^t  is  located  at  its 
corresponding  nominal  point  described  abo*  .. 

The  exponent  of  the  joint  Gaussian  distribution  Is  a  quad¬ 
ratic  form  and  represents  an  ellipsoid  in  the  h»  V»  yapskce.  The 
orientation  and  semi-major  axes  for  this  ellipsoid  are  illuminated 
by  a  coordinate  rotation  of  this  quadratic  form  into  its  normal  co¬ 
ordinates.  such  that  all  cross-product  terms  of  the  exponent  dir 
appear.  In  terms  of  this  new  set  of  coordinates,  the  variables  h'» 

V',  Y*  are  mutually  and  pairwise  independent.  Although  not  generally 
true,  it  is  true  for  the  Gaussian  distribution  that  if  the  correlations 
between  variable#  are  zero,  the  variables  are  independent.  Hence, 
the  distribution  in  normal  coordinates  is  a  joint  function  of  three 
independent  variables. 

Since  they  are  expressed  in  terms  cf  zero  means  and  nor¬ 
malized  relative  to  their  respective  variances,  the  three  variables 
thus  reduce  to  three  independent  normal  variates  of  zero  mean  and 
unity  variance.  It  is  an  established  fact  that  the  distribution  of  the 
sum  of  the  squares  of  n  independent  normal  variables  of  zero  mean 
and  unit  variance  is  the  chi-squared  distribution  with  n  degrees  of 
freedom.  Hence,  after  a  rotation  of  the  quadratic  form  to  normal 


where 

l-^lmn 

ym*y„ 
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coordinates  the  joint  distribution  is  reduced  to  a  chi- squared  dis* 
trlbution  of  3  d.f. 

Let  Q  be  the  symmetrical  matrix  of  the  original  quadratic 
form.  Then  Q  can  be  reduced  to  diagonal  form  by  some  normalized 
orthogonal  matrix 

(29)  CtX  =Xk 


where  K  is  the  diagonal  matrix  of  eigen  values  for  Q.  It  remains 
to  solve  for  and  Q. 

The  solution  of  equation  (2))  leads  to  a  set  of  simultaneous 
equations 

<?11  "^1'^  11  + ^21  ■*■•••  ^In^nl  =  ^ 


(30) 


^21^12+  <^21  ■’'2>^22  +  -  •  •  »2n^n2 


^  *nl^ln  +  »n2^2n  +  -  •  • 


nn  n'  nn 


where  a.^  =  element  of  Q 
f  =  element  of^^ 

k.  =  element  of  K 

This  set  of  equations  has  a  non'-trivial  solution  only  if 


(31) 


ail  -k 

.ava. 

12  In 

^21 

*22  ■  ^  ...  “2n 

^nZ 

-  *nn  - 

which  leads  to  an  n^  order  plynomxal  in  k.  In  matrix  notation  this  reads 

[q  -  k|  =  0. 


26 


Thus,  from  equation  (27)  the  quadratic  form  is 
(32)  lAJhh  +  l-^VV  +l^lvY 

-  zlAl^yhV- 2|A|^^  W  -2lAj^^Vh=)c 

where  the  constant  k  is  established  in  order  to  make  this  ellipsoidal 
volume  a  99  percent  (3^)  confidence  interval. 

The  direction  numbers  for  this  quadratic  form  are  established 
by  means-of  the  following  set  of  equations,  i  b  1.  2.  3. 


where  the  k^  are  the  eigen  values  for  the  quadratic  form  and  are 
computed  from  the  third  degree  polynomial. 


Nhh-. 

|A| 

hV 

I^'hy 

1^1  hV 

1A|  yy  -  k. 

II 

lAlhv 

Hence  for  each  k..  a  triple  ratio  ^y/m^/n^  6i  direction  numbers  is 
obtained,  where  each  such  act  defines  the  direction. of  one  of  the 
normal  coordinates. 
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From  V  chi-^squared  distribution  table*  it  is  found  that  k  z  11.341 
for  3  d.f.  and  the  quadratic  surface  .tn  nominal  coordinates  (hj#  Vj^,  Yj^)  is 

(35)  kjh^  +•'2^.^  +J'3V,f  =  ”*341 

or 


(36) 


Thus*  the  99  percent  confidence  ellipsoid  is  defined,  it  has  a 
canonic  form  relative  to  the  coordinate  system  defined  by  equation  (37): 

(“ 

(37)  n>, 

\ 

and  the  semi-major  axes  in  terms  of  3<r  values  are 


(-^y  .(M-  ,  Qd 

-11.341  +  ~1?:341'  •  *-  iiTTir 

~9Vc,  9  k,  9  k, 


=  1 


(38) 

<3'hn>^ 

11.341 

-  Ok, 

(39) 

(3.v„)2 

11.341 

(40) 

11.341 

=  -^E7“ 

As  an  illustration  of  the  above  theory,  it  will  be  applied  to  the  two- 
dimensional  case  where  y  =  0; 
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2  2 


(41) 


p(h.V)  = 


where  =  nij^  =  0 


(42)  A  = 


••hV  Vv 


Pvh  “’n'^V 


(43) 


I'N 


2  2  2  2  2  _2_2,,  ,2, 

“^V  '  •’hv  “^h  ''v  ■  h  “^V  '*  "  Phy  ' 


O  a  Quadratic  form 


(44)  Q  =  — 2— y 


2^h‘<rv‘(l-PhV> 

N 


lAl^vhV  +  lAlvH  Vlx  +  lAlvvV^J 


hh  ’ 

hv  =  -Pvh  Vv 


(45) 


lAl 

|a! 


V.  =  "PhyVv 

n 


VV  “  ''h 


(46) 


'  2,  T—  <^y  -  2p  +  “-h  V  » 


2(1  -  p‘)  <r‘  0-, 


h  “V 
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Solving  lor  the  eigen  values,  in  terms  of  the  determinant  for 
the  matrix  o£  cofactors, 

2 


(47) 


-k 
-p  Vv 


-p 

2  , 


=  0 


(48) 


.2^  2  ^  [Ti  2.2.  . _ ^2  2 

(^h  +  "'v  *  /y  ■  "^v'  P 

- 2 - 


Now  solve  for  the  normal  directions  of  the  ellipse: 


(50) 


I  (Ty  -  k)  f 
,_-P 


(51) 


'(■Ty  -  k)  f 
P  ThO-yl 


+ 


p  =  0 

i<r^  -  k)m  =  0 


for  k  =  kj,  k.^ 


4p  o-j^CTym 
+(<rj^  -  k)m 


(52) 


i  /m  = 


+p  0-.  flr. 


h''V 


Oy  - 


<^h 

p  Vv 


Either  relation  will  lead  to  the  same  result.  The  mrection 
cosines  for  the  normal  coordinates  are  therefore  the  elements  of  the 
following  transformation  matrix. 
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v/here'h  and  V  are  components  of  a  vector  relative  to  the  original 
coordinates,  and  h'  and  V*  are  the  components  of  the  same  vector 
relative  to  the  normal  coordinates. 


Define : 


(54) 


^  2 

TiT 


“‘2 

rr 


‘1  _ 
^  - 


P  O-hO-. 


h"V 


/Tz  ,  ,2  2  2  2 

•V  ■  *'1^  P  “’V 


''y  -  I 


-  k,)2  +  0- 


-h  ->'2 


V‘^h  +  P 


2  _  2 
V 


po-hO-y 


‘^h  ""v 
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Hence 

(55) 


ta„e  = 

TJ  P'h”'v 


for  a  counterclockwiac  rotation  from  the  V  axes  (see  Figure  8). 


2  2 

2  . 

■^h  ■  *'2  -  r  ■  T 


2,2  .  2  2  2 
,)  +4p  <r^  fl-y 


where 


,  ,  /TZ  2.2  ^  .2  Z~Z 


Hence 


(56) 


l.et 


therefore 


tan  6=0 


tan  26  : 


-  I  '"h  1:!^  +  Sldtl  ,  1 

h^v  j^o-h^v  'V^^vf.Ty^  j 


2  tan  0 

- ^2 - 

I  -  tan^  0 


<'h  -  -V ) 


■  Vv 

tan  e  =  -  ^  ^ 


and 


tan  20  = 


-2(a+fyaZ  +  l 


-(2)(a^  +  aVa^  +  1) 
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thus 


<57  ) 


Define: 


-1 


Therefore  ^ 


( 


or 


2 

V 


Pvh 


since  0^  is  a  normalized 
orthogonal  matrix 


Hence 

(58,  Q.  =  -^  +k2vfj 


For  2  d.f.  we  find  rom  the  chi- squared  table  that  the  ex¬ 
ponent  k  should  be  9>937  in  order  *hat  the  dispersion  ellipse  include 
99  percent  of  all  dispersions  from  the  nominal.point  at  main  retro 
burnout.  Thus# 


33 


(59) 


9.937 


k,  hf /(9)  k2  VjV9 

iPSiimr  *  T;9377(5r 

Divide  both  numerator  and  denominator  by  9  in  both  term* 
to  force  the  aemUmajor  axes  of  the  dispereion  ellipse  to  be  Bcr  values. 
Hence,  the  semi-major  axes  about  the  nominal  point  are; 

(60)  (3(rvj)  =  — 
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2.  MAXIMUM  VERNIER  THRUST  SURFACE 

The  maximum  vernier  thrust  surface  is  defined  as  the  sur¬ 
face  of  points  in  the-h.  V.  v  apace*  from  which  the  vehicle  can 
come  to  a  zero'altitude  zero  velocity  landing*  if  retroed  at  the  znaxi- 

mum  vernier  thrust.  The  equations  of  motion  for  this  phase  of  the  T 

trajectory  are  (assuming  constant  mass): 


(58) 

&  1) 

a 

(59) 

fix  *  “  8  b 

3 

(60) 

Vy  Sy  dt  -  Vg  sin  yg 

S 

es 

(61) 

^  +  V® 

pc 

U) 

(62) 

Vx  -  -X  ’'2 

pt 

el 

(63) 

y  Vydc 

/>t 

t6 

(64) 

X  -  j  dc 

(6 

a^*i=x,y  Acceleration*  ft/sec^ 
g  Lunar  gravity,  ft/sec^ 

b  Maximum  vernier  thrust/lunar  weight  ratio 

E< 

sc 

Vj,  isx,  y  Velocity,  It/sec 

(h 

::*y  Displacement,  feet 

ti( 

pi 

This  set  of  equations  is  solved  simultaneously  until  0. 

At  this  time*  h2  =  **/•  In  this  way  a  number  of  discrete  points  on 

the  surface  are  determined.  Then*  by  means  of  a  curve  fitting  pro-  (6 

cess*  a  smooth  curve  can  be  found  to  satisfy  all  points  found. 

For  the  data  run  to  date*  a  surface  that  approximates  all  (6 

data  points  to  within  3  percent  is 
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The  time  for  descent  is  found  in  a  similar  fashion  and  is 


3.  SURFACE  OF  NOMINAL  MAIN  BURNOUT  CONDITIONS 

Having  defined  the  distribution  ellipse.  E»  and  the. surface, 
in  the  hj  Yj  space,  a  surface  of  nominal  points  can  be 
established.  This  surface  is  a  locus  of  main  retro  burnout  points 
under  ideal  conditions,  i.e..  zero  error  perturbations. 

This  surface  is  established  by  finding  the  centers  of  the 
ellipsoids  which  are  tangent  to  the  surface. 

VpS/g 

(I  +  1/810  yg)-i 

''w.  -  •'n)®  *  Aw  <''a  -  (7a  - 

(6V;  .  3  A^v  ■  ”»)  ♦  2*h7  '  ’’nXra  -  7n) 

*  a  Avy  (Vg  -  v„)(7a  -  7„)  .  i.2£  .  c 


Equations  (66)  and  (67)  represent  the  maximum  vernier  thrust 
surface  and  the  elliptical  surface  in  terms  of  a  nominal  center. 

(h  ,  V  .  Y  )*  thc&r  point  of  tangcncy,  (h^*  Y^)*  ^  addi¬ 

tion.  the  slopes  of  the  two  surfaces,  in  the  h-V  plane  and  h-Y 
planes  must  also  be  equal.  These  latter  relations  arc  described 
by  equations  (68)  and  (69). 

2hp  Ayy  (Vp  -  V„)  .  A^y  (h  p-  ilp)  *  A  ( /g  -  7„ ) 


Ahh  *  AhV  (Vg  -  v„)  +  A^^  (rg  -  7„) 


J7-5  ^  su? 


(rp  -  /n)  *  AhY  (!■.  -  l>n)  ^  Ay.  (Vp  -  Vn) 
(hg  -  I'n)  *■  Ahv  (Vg  -  V„)  ♦  Aj,y  (rg  -  7n' 
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Equations  (66).  (67).  (68).  and. (69)  are  solved  simultaneously  by 

assuming  values  for  and  Y2*  solving  for  b^.  V^.  and  y^ 

under  the  condition  that  h  >  h,. 

n 

The  symbols  employed  here  have  been  defined  previously. 
The  surface  of  nominal  points  can  then  be  established  by  means  of  a 
curve  fitting  procedure*  or  the  data  points  can  merely  be  stored  as 
such  for  future  use.  This  surface  (or  points)  shall  hereafter  be  re¬ 
ferred  to  as  the  surface. 

4.  SURFACES  OF  CONSTANT  VERNIER  FUEL 


The  surface  of  constant  vernier  fuel  is  established  by  inte¬ 
grating  the  equations  of  motion  during  the  drift  period  (minimum 
vernier  thruat)  backwards  (negative  time)  using. the  points  of  S^  as 
initial  conditions. 


(70) 


<=-'"2*2 


(71) 

(72) 


t|  e  time  for  minimum  vernier  thrust  phase*  seconds 
t^  =  time  for  maximum  vernier  thrust  phase* 
seconds  (see  equation  (65)) 
m^  =  minimum  thrust  mass  flow  rate*  slugs/sec 
(a  constant) 

m2  =  maximum  thrust  mass  flow  rate*  slugs/sec 
(a  constant) 

C  =  total  vernier  fuel  mass*  slugs 

r^i  V 

Vy  ■■  j  8  (a  y  +  l)  dt  -  Vg  sin  yg 

8  a  “  dt  +  Vg  cos  72 


(73) 

(74) 


V  +  Vy2 

r‘*  .  . 

h  •*  ;  -  Vw  dt  rto 
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(75) 


where  is  defined  by  equation  (64). 

For  every  va'lue  of  hj.,  and  y>  a  is  defined,  (see 
equation  (65)).  Hence,  every  value  C  establishes  a  tj,  a  minimum 
vernier  drift  time.  Integrating  equations  (70)  through  (75)  yields 
points  (hj,  V^,  Yj),  for  various  values  of  C,  By  curve  fitting  these 
points  $£  is  established,  with  C  asa  parameter. 


(76) 
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la  -  a)  ^  4  a  V]  -  V-^y 


Wp  Si  vernier  fuel  weight 

=  average  vehicle  weight  during  vernier  descent 


where 

(77) 

b'  =  b  -  0.57  (1  -  sin  Yf) 

(78) 

V'  =  t  1  -  0.38  (1  -  sin  Yf)  3 

5.  INTEGRATION  OF  THE  COMPUTER  PROGRAM 

The  IBM  programs  for  the  various  surfaces: 

£  s  dispersion  ellipsoids 
Sj^  s  surface  of  nominal  points 

=  surface  of  maximum  vernier  fuel 
=  surface  of  constant  vernier  fuel  , 

are  integration  into  a  master  program.  This  master  program 
accepts  any  set  of  perturbations  and  initial  conditions  as  inputs 
and  establishes  the  various  surfaces  defined  above. 

The  total  fuel  is  then  optimized  as  a  function  of  main  retro 
terminal  velocity  in  the  following  manner: 

1)  A  set  of  main  retro  terminal  conditions  establishes  a 
dispersion  ellipsoid  center  and  surface. 

2)  A  computer  search  is  then  made  for  the  point  ofvtan- 
gency  of  this  ellipsoid  and  the  surface.  Ti\is  is  done 
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by  finding  the  intersections  of  £  and  for  various 
values  of  and  C.  When  an  surface  is  found  (identi¬ 
fied  by  C)  such  that  its  intersection  with  the  ellipsoid 
is  one  point,  the  search  ends.  This  point  of  i.angency 
represents  the  maximum  required  amount  of  vernier 
fuel  rec  a*rcd  corresponding  to  the  chosen  nominal 
point. 

3)  The  totax  engine  weight  is  then  computed  as 


(79) 


k  +  K 
mm  V 


C  = 


a  ratio  of  total  main  retro  engine  weight  to  its 
propellant  weight 

Ky  a  ratio  of  total  vernier  engine  weight  to  its 
propellant  weight 
s  total  engine  weight 

4)  Curves  of  k^M _ and  K,,  C  versus  nominal  main  retro 

m  m  V 

terminal  velocity  V.J.,  are  generated  by  repeated  com¬ 
puter  runs  for  a  range  of  desirable  values  for  V.^.  The 
sum  of  the  two  engine  weights  is  characterized  by  a 
minimum  at  some  nominal  main  retro  terminal  velocity. 
This  value  of  terminal  velocity  represents  the  optimum 
in  the  trade-off  of  energy  removal  during  the  main  retro 
and  verniar  phases  as  discussed  at  the  end  of  Part  1. 

The  above  program  is  integrated  so  that  a  given  set  of 
perturbations  and  initial  conditions  will  result  in  an 
engine  system  design  of  least  weight. 
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CONCLUSION 


Zn  summary*  broad  considerations  of  simpUcity-and.seliabslity, 
and  hardware  trade-offs  involving  accuracy  and  weight*  have  led  to 
the  choice  of  a  two-enginc  relxo  system  for  the  slowdown  from  lunar 
approach  speed  to  a  soft  landing  speed.  Within  the  system  limitations 
imposed  by  the  practical  propulsion  and  sensing  techniques  existing 
today*  an  optimization  procedure  has  divided  the  energy  removal 
between  the  two  engines  such  that  total  retro  system  weight  is  least. 

A  digital  computer  program  for  the  optimization  in  the  case  of  off- 
vertical  trajectories  has  also  been  detailed. 

Should  the  available  engine  and  sensor  techniques  change*  and 
they  will*  the  answers  shown  and  implied  here  will  also  change.  But 
it  seems  that  the  broad  design  approach  and  the  optimization  procedure 
will  remain  ''alid  and  that  their  use  with  different  inputs  will  produce 
new*  correct*  different  results. 
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FIGURE  4-  VERNIER  ENGINE  DESCENT 
AND  VERNIER  FUEL 


NOMINAL  MAIN  RETRO  BURNOUT  VELOCITY 


WEIGHT 


FIGURE  5  CHANGE  IN  FUEL  AND 

system  weights 


